Critical behavior at Mott-Anderson transition: a TMT-DMFT perspective 
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We present a detailed analysis of the critical behavior close to the Mott-Anderson transition. 
Our findings are based on a combination of numerical and analytical results obtained within the 
framework of Typical-Medium Theory (TMT-DMFT) - the simplest extension of dynamical mean 
field theory (DMFT) capable of incorporating Anderson localization effects. By making use of 
previous scaling studies of Anderson impurity models close to the metal-insulator transition, we 
solve this problem analytically and reveal the dependence of the critical behavior on the particle-hole 
symmetry. Our main result is that, for sufficiently strong disorder, the Mott-Anderson transition 
is characterized by a precisely defined two-fluid behavior, in which only a fraction of the electrons 
undergo a "site selective" Mott localization; the rest become Anderson-localized quasiparticles. 
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Many strongly correlated materials find themselves 
close to Mott localization [l| - a process through which 
all the valence electrons within a narrow band turn into 
localized magnetic moments. In real systems, disorder 
introduced by doping or impurities often cannot be ne- 
glected, as it provides an alternative fundamental mech- 
anism for suppressing metalicity, through the process of 
Anderson localization [2|. The effects of weak interac- 
tions in this regime have been studied using perturbative 
methods [3|, but these approaches cannot describe the 
strong correlation effects associated with incipient mag- 
netism and Mott localization. 

Which of these two routes to localization - Anderson 
or Mott - dominate? In most cases, simple estimates 
show that both effects play a comparable role and both 
need to be taken into account. Most existing theories are 
not able to combine these two fundamental processes in 
the same framework, and this conceptual difficulty has 
provided the essential pitfall in our understanding of the 
metal-insulator transition (MIT). 

At the moment, the most successful theory for the 
Mott transition is based on dynamical mean field the- 
ory (DMFT) [4[ ideas. By replacing the environment 
of each site by its average value, the original version of 
this theory proved unable to describe the spatial fluctu- 
ation effects associated with the approach to the Ander- 
son transition. Very recent work [5j], however, identified 
the conceptually simplest extension of DMFT capable to 
overcome these shortcomings - the Typical- Medium The- 
ory (TMT-DMFT). In the non-interacting limit, this the- 
ory provides a reasonable picture of the Anderson transi- 
tion, as established by quantitative comparison [5| with 
exact (numerical) results. 

The TMT-DMFT method was first applied to the dis- 
ordered Hubbard model by Byczuk et al. @, who ob- 



tained the phase diagram for this problem from the nu- 
merical solution using the Numerical Renormalization 
Group (NRG) method for the impurity solver. However, 
the physical nature of the phases and of the phase tran- 
sition was not investigated in that numerical study. 

The task of elucidating the physical mechanism and 
the precise form of the Mott-Anderson critical point 
within the TMT-DMFT description is the main subject 
of this Letter. By making use of previous scaling stud- 
ies pj of Anderson impurity models close to the MIT, 
we present a detailed analytic solution for this problem, 
which emphasizes the dependence of the system proper- 
ties on its particle-hole symmetry. Our main finding is 
that, for sufficiently strong disorder, the physical mecha- 
nism behind the Mott-Anderson transition is the forma- 
tion of two fluids, a behavior that is surprisingly reminis- 
cent of the phenomenology proposed for doped semicon- 
ductors [|| . Here, only a fraction of the electrons (sites) 
undergo Mott localization; the rest can be described as 
Anderson-localized quasiparticles. Thus, in our picture 
the Mott-Anderson transition can be seen as reminiscent 
of the "orbitally selective" Mott localization [9( ; precisely, 
here we have a "site selective" Mott transition, since it 
emerges in a spatially resolved fashion. 

TMT-DMFT and order parameters - We consider a 
half- filled Hubbard model [J] with random site energies, 
as given by the Hamiltonian 



H = -V 



<ij>a 



C ia C 3 a 



£in i<7 + U^n^nii, (1) 



Here, c\ a (cj CT ) creates (destroys) a conduction electron 
with spin a on site i, rii a = c\ a Ci a , V is the hopping 
amplitude, and U is the on-site repulsion. The random 
on-site energies Si follow a distribution P(s), which is 
assumed to be uniform and have width W. 
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TMT-DMFT @, @] maps the lattice problem onto 
an ensemble of single-impurity problems, correspond- 
ing to sites with different values of the local energy 
£j, each being embedded in a typical effective medium 
which is self-consistently calculated. In contrast to stan- 
dard DMFT [nj, TMT-DMFT determines this effective 
medium by replacing the spectrum of the environment 
("cavity") for each site by its typical value, which is de- 
termined by the process of geometric averaging. For a 
simple semi-circular model density of states, the cor- 
responding bath function is given by @, @] A(w) = 
V 2 Gtyp(cL>), with G typ (uj) = j^duj 1 p typ {uj')/{uj-uj') be- 
ing the Hilbert transform of the geometrically-averaged 
(typical) local density of states (LDOS) p t y P (ui) = 
exp{J deP(e) In p(u), e)}. Given the bath function A(u>), 
one first needs to solve the local impurity models and 
compute the local spectra p(ui,e) = — n^ 1 ImG(w, e), 
and the self-consistency loop is then closed by the the 
geometric averaging procedure. 

To qualitatively understand the nature of the critical 
behavior, it is useful to concentrate on the low-energy 
form for the local Green's functions, which can be speci- 
fied in terms of two Fermi liquid parameters as 



G(w,e,) 



ZiA(uy 



(2) 



where Zi is the local quasi-particle (QP) weight and ii is 
the renormalized site energy • The parameters Zj and 
£j can be obtained using any quantum impurity solver, 
but to gain analytical insight here we focus on the varia- 
tional calculation provided by the "four-boson" technique 
(SB4) of Kotliar and Ruckenstein ll|, which is known to 
be quantitatively accurate at T = 0. We should stress, 
though, that most of our analytical results rely only on 
Fermi liquid theorems constraining the qualitative be- 
havior at low energy, and thus do not suffer from possible 
limitations of the SB4 method. 

Within this formulation, the metal is identified by 
nonzero QP weights Zi on all sites and, in addition, 
a nonzero value for both the typical and the average 
[p av (ui) — J deP(e)p(ui,e)] LDOS. Mott localization (i.e. 
local moment formation) is signaled by Zi — ► [l0(, 
while Anderson localization corresponds to Z ; / and 
p av ^ 0, but p typ = @, |||. While Ref. [|| concentrated 
on ptyp and p av , we find it useful to simultaneously ex- 
amine the QP weights Zi, in order to provide a complete 
and precise description of the critical behavior. 

Phase diagram - Using our SB4 method, the TMT- 
DMFT equations can be numerically solved to very high 
accuracy, allowing very precise characterization of the 
critical behavior. In presenting all numerical results we 
use units such that the bandwidth B = 4 V = 1. Fig. la 
shows the resulting T — phase diagram at half filling, 
which generally [13] agrees with that of Ref. [f| . By con- 
centrating first on the critical behavior of the QP weights 
Zi, we are able to clearly and precisely distinguish the 



»2 





£ 




-1.6 -0.8 0.0 


0.8 1 


(a) 


V 1 


T — T — 1 


Mott-Anderson insulator 






(gapless) 


,\ 


l| 


© &■ q / 

■ y \ / 


l\ 

f\ 


ft 

mi 
/I' 
Jt, 


1 d> 

{ * 

Correlated metal / 


(c) 




t 

2 Mott insulator 
1 (hard gap) 








^ , - 






i ■ A 1 


.-)—.■-!■ 


T -I-'.'- 



-0.6 -0.3 0.0 0.3 0. 



FIG. 1: (Color online) (a) T = phase diagram for the disor- 
dered half filled Hubbard model, obtained from the numerical 
SB4 solution of TMT-DMFT. Panels (b) and (c) show the evo- 
lution of the quasi-particle weight Z(si) in the critical region. 
Behavior at (b) the Mott-Anderson transition (W > U) is 
illustrated by increasing disorder W = 2.5, 2.6, 2.7, 2.8, 2.83 
(from the black curve to the brown one), for fixed U = 1.25; 
and at (c) the Mott-like transition (IT < U) by increasing 
the interaction interaction U = 1.5, 1.6, 1.7, 1.8 and 1.86, at 
fixed disorder W = 1.0. 



metal from the insulator. We find that at least some of 
the Zi vanish all along the phase boundary. By taking 
a closer look, however, we can distinguish two types of 
critical behavior, as follows. 

Mott-Anderson vs. Mott-like transition - For suffi- 
ciently strong disorder (W > U), the Mott-Anderson 
transition proves qualitatively different than the clean 
Mott transition, as seen by examining the critical behav- 
ior of the QP weights Zi = Z{ei) (Fig. lb). Here Z t -> 
only for < |ej| < U/2 , indicating that only a fraction of 
the electrons turn into localized magnetic moments. The 
rest show Z^ — * 1 and undergo Anderson localization (see 
below) . Physically, this regime corresponds to a spatially 
inhomogeneous system, with Mott fluid droplets inter- 
laced with regions containing Anderson-localized quasi- 
particles. In contrast, for weaker disorder (W < U) the 
transition retains the conventional Mott character. In 
this regime Z.- L — > on all sites (Fig. lc), corresponding 
to Mott localization of all electrons. 

Wavefunction localization - To more precisely charac- 
terize the spatial fluctuations of the quasiparticle wave- 
functions, we compare the behavior of the typical {ptyp) 
and the average (p av ) LDOS. The approach to the Mott- 
Anderson transition (W > U) is illustrated by increasing 
disorder W for fixed U = 1.25 (Fig. 2 - top panels). 
Only those states within a narrow energy range (u < t, 
see also Fig. 4) around the band center (the Fermi en- 
ergy) remain spatially delocalized (pt VP ~ Pav), due to 
strong disorder screening 0, within the Mott fluid 
(sites showing Zi — > at the transition). The electronic 
states away from the band center (i.e. in the band tails) 
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FIG. 2: (Color online) Frequency dependence of p typ (full 
line) and p av (dashed line) in the critical region. Results in 
top panels illustrate the approach to the Mott- Anderson tran- 
sition (W > U) at U = 1.25; the bottom panels correspond 
to the Mott-like transition (W < U) at W = 1.0. 



FIG. 3: (Color online) Typical and average values of p(0) as 
the metal-insulator transition is approached for (a) U — 1.25 
and (b) W = 1.0. The insets show p(0) as a function of e for 
(a) W = 2.5, 2.6, 2.7 and 2.83 (from the black curve to the 
blue one) and (b) U = 1.5, 1.6, 1.7 and 1.86. 



quickly get Anderson-localized, displaying large spatial 
fluctuations of the wavefunction amplitudes [5| and hav- 
ing ptyp < Pav 

The spectral weight of the delocalized states (states in 
the range oj < t) decreases with disorder and vanishes 
at the transition, indicating the Mott localization of this 
fraction of electrons. At this critical point, the crossover 
scale t also vanishes. In contrast, the height ptyp(0) re- 
mains finite at the transition, albeit at a reduced In- 
dependent value, as compared to the clean limit. More 
precise evolution of pt yp (0) is shown in Fig. 3a, demon- 
strating its critical jump. 

Behavior at the Mott-like transition (W < U) is 
dramatically different (Fig. 2 - bottom panel). Here 
Ptyp ~ Pav over the entire QP band, indicating the 
absence of Anderson localization. It proves essentially 
identical as that established for the disordered Hub- 
bard model within standard DMFT [l(| , reflecting strong 
correlation-enhanced screening of disorder Q, where 
both p av (j-o = 0) and p typ (aj = 0) approach the bare 
(W = 0) value (see also Fig. 3b). Similar results were 
found in Ref. 0], but an explanation was not provided. 

The corresponding pinning @, EH for p(uj — 0, e) is 
shown in the insets of Fig. 3, both for the Mott- Anderson 
and the Mott-like transition. In the Mott-Anderson 
case, this mechanism applies only within the Mott fluid 
(|e| < U/2), while within the Anderson fluid (|e| > U/2) 
it assumes smaller values, explaining the reduction of 
Pty P (0) in this case. We suggest that this spatial distribu- 
tion of the DOS at the Fermi energy (each e corresponds 
to a different position in the lattice) could be probed by 
scanning tunneling microscopy experiments. 

Analytical solution - Within our SB4 approach, the 
TMT-DMFT order-parameter function p ty p(ui) satisfies 



the following self-consistency condition 

Pty P (uj) = exp J deP(e) {\n[V 2 Z 2 (e)p typ (uj)} 

- ln[(u> - e{e) - V 2 Z(e) Re G ty p{uo) f 

+ (nV 2 Z(e)p t yp(u;)) 2 }}. (3) 

While the solution of this equation is in general difficult, 
it simplifies in the critical region, where the QP param- 
eter functions Z{e) and s{e) assume scaling forms which 
we carefully studied in previous work Q. This simplifi- 
cation allows, in principle, to obtain a closed solution for 
all quantities. In particular, the crossover scale t, which 
defines the p t yp{w) mobility edge (see Fig. 4 and Ref. 0]), 
is determined by setting ptyp (u = t) = 0. 

Using this approach we obtain that, in the case of 
Mott-like transition (W < U), the critical behavior of all 
quantities reduces to that found in standard DMFT [l(| , 
including t ~ U C (W) — U (in agreement with the nu- 
merical results of Fig. 4b), perfect screening of site ran- 
domness 0, and the approach of p av (^ = 0) and 
Ptyp(tu — 0) to the clean value. The precise form of the 
critical behavior for the crossover scale t is more com- 
plicated for the Mott-Anderson transition (W > U) (as 
confirmed by our numerical results in Fig. 4a), and this 
will not be discussed here. 

Instead, we focus on elucidating the origin of the puz- 
zling behavior of p c = pt yp {u — 0), which is known || 
to vanish linearly p c ~ (W c — W) for U — 0, but which 
we numerically find to display a jump (i.e. a finite value) 
at criticality, as soon as interactions are turned on. For 
lu = our self-consistency condition reduces [l3} to 

h p ^ m>?fv*U)Pi =0 ' (4) 

which further simplifies as we approach the critical point. 
Here, the QP parameters Z(e) — > and e(e) ~ Z 2 (e) -€i 
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FIG. 4: (Color online) Frequency dependence of the typical 
DOS very close to the metal-insulator transition for (a) the 
Mott-Anderson transition (W > U) at U = 1.25 and (b) 
the Mott-like transition (W < U) at W = 1.0. The insets 
show how, in both cases, the pt yp {w) bandwidth t — > at the 
transitions. 



Z(e) for the Mott fluid (|e| < U/2), while Z(e) — > 1 and 
\e(e)\ — > \e - U/2\ for the Anderson fluid (|e| > U/2), 
and we can write 



= 



f U/2 1 

/ deP(e)ln 

Jo (nV t 



{nVpcY 

(W~U)/2 

deP (s) H(e/V) 2 + (wVp c 



(5) 



This expression becomes even simpler in the U << W 
limit, giving 

^-ln-?—+a-bVp c + O{pt)=0, (6) 

where a(W, U) = {l-U/W){l-hx[{W-U)/2V]} and b = 
. This result reproduces the known result [B[ p c ~ 
(W c — W) at U = 0, but dramatically different behavior is 
found as soon as U > 0. Here, a non-analytic (singular) 
contribution emerges from the Mott fluid (|e| < U/2), 
which assures that p c must remain finite at the critical 
point, consistent with our numerical results (see Fig. 3). 
Note that the second term in Eq. ([5|), coming from the 
Anderson fluid (|e| > U/2), vanishes in the case of a 
Mott-like transition (U > W), and our result reproduces 
the standard condition np c V = 1 [Toj | . which corresponds 
to the clean limit. 

A further glimpse on how the condition Ttp c V = 1 is 
gradually violated as we cross on the Mott-Anderson side 
is provided by solving Eq. ([5]) for U ^ W limit, giving 



1 



24 \V 



u_ 

w 



(7) 



again consistent with our numerical solution [14j. 

But what is the physical origin of the jump in p c ? 
To see it, note that the singular form of the first term 



in Eq. {5} comes from the Kondo pinning [10| e(e) ~ 
Z 2 (e) < Z(e) within the Mott fluid. This behavior re- 
flects the particle-hole symmetry of our (geometrically 
averaged) ptyp(y> = 0) bath function, which neglects site- 
to-site cavity fluctuations present, for example, in more 
accurate statDMFT theories [l5j . Indeed, in absence of 
particle-hole symmetry, one expects [Io[ ~ Z(e), and 
the resulting e-dependence should cut-off the log singu- 
larity responsible for the jump in p c . This observation 
provides a direct path to further refine the TMT-DMFT 
approach, reconcilin g th e present results with previous 
statDMFT findings [l5(. As a next step, one should 
apply the TMT ideas to appropriately chosen effective 
models [l6(, in order to eliminate those features reflect- 
ing the unrealistic particle-hole symmetry built in the 
current theory. We emphasize that the two-fluid picture 
is a consequence of only a fraction of the sites showing 
Z —> and is not dependent on cither particle-hole sym- 
metry or the consequent jump in the DOS. 

Conclusions - This Letter explores the TMT-DMFT 
critical region of the Mott-Anderson transition. We show 
how key insight can be obtained by focusing on the evo- 
lution of the local quasiparticle weights Zi as a second 
order parameter describing tendency to Mott localiza- 
tion, in addition to the Anderson-like TMT order pa- 
rameter ptyp- This analysis reveals the fundamental two- 
fluid character of the Mott-Anderson transition, consis- 
tent with the phenomenology proposed for doped semi- 
conductors [S\. Physically, it describes spatially inhomo- 
geneous situations, where the Fermi liquid quasiparticles 
are destroyed only in certain regions - the Mott droplets 
- but remain coherent elsewhere. Understanding the de- 
tails of such "site selective" Mott transitions should be 
viewed as an indispensable first step in solving the long- 
standing problem of metal-insulator transitions in disor- 
dered correlated systems. 
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